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ABSTRACT 



Aims. We present an adaptation of the standard scenario of disk-galaxy formation to the concordant ACDM cosmology aimed to derive 
analytical expressions for the scale length and rotation speed of present-day disks that form within four different, cosmologically 
motivated protogalactic dark matter halo-density profiles. 

Methods. We invoke a standard galaxy-formation model that includes virial equilibrium of spherical dark halos, specific angular 
momentum conservation during gas cooling, and adiabatic halo response to the gas inflow. The mean mass-fraction and mass-to-light 
ratio of the central stellar disk are treated as free parameters whose values are tuned to match the zero points of the observed size- 
luminosity and circular speed-luminosity relations of galaxies. 

Results. We supply analytical formulas for the characteristic size and rotation speed of disks built inside Einasto r''*, Hemquist, 
Burkert, and Navarro-Frenk- White dark matter halos. These expressions match simultaneously the observed zero points and slopes 
of the different correlations that can be built in the RVL space of disk galaxies from plausible values of the galaxy- and star- formation 
efficiencies. 

Key words, dark matter - Galaxies: formation - Galaxies: fundamental parameters - Galaxies: spiral - Galaxies: structure 



in 

u 

o 

-(— > 

L^ 1. Introduction 



l> 

p 

o 
o 



X 



In the current hierarchical galaxy-formation paradigm disk- 
galaxies are born out of the hot gas-atmospheres associated with 

[ the potential well of virialized cold dark matter (CDM) halos. It 
is assumed that baryons have initially both the same density pro- 
file and specific angular momentum distribution as DM - the lat- 
ter achieved, for instance, through tidal inter actions vyith ne igh- 
boring objects in the precollapse phase (e.g.l PeeblesLll969l) . As 

' the gas radiates its energy it cools and starts to fall towards the 
center of the DM halo maintaining its specific angular momen- 
tum, where it settles into a rotationally supported disk. The as- 
sembly of a concentration of cold baryons at the bottom of the 
gravitational potential well on timescales longer than the free- 
fall time produces the adiabatic contraction of the dark halcQ. In 
this standard picture, the internal properties of disk galaxies are 

< expected to be largely dictated by those of their host halos, and 
through the latter, by those of the background cosmology too. 

Theoretical predictions for the distribution of disk galax- 
ies in the space of disk scalelength (or size), fiducial (usually, 
maximum or asymptotic) rotational speed, and luminosity (or 
mass) based, partially or totall y, on the scenari o just outl ined are 
abundant in the literature (e.g.jMo et al.|[T998l hereafter lMMWt 
iPiza^no et al J 1200 5: Dutton et al. '2007). They are widely used 
in semi-analytic cosmological models, pre-prepared numerical 



' In modern literature, the mode and amount of ha l o contraction 
are ac tually a matter of debate (e.g. iDutton et al.L 1 20071 : iTissera et al.L 
120091) . The outcome, however, remains unchanged: the properties of 
disk galaxies are linked to those of their host halos. 



simulations of galaxy groups and clusters, and studies of disk- 
galaxy scaling relations. 

While nowadays there are extensive and comprehensive in- 
vestigations of the correlations between disk-galaxy properties 
that deal with the scatter and covariances of the variables and 
allow for different modes of halo contraction (e.g. lDuttonet all 
i2007), it is not always feasible to implement such sophisticated 
treatments whenever one needs to estimate the scaling of the ba- 
sic structural and kinematic parameters of galaxies. The simplest 
alternative is the use of scaling laws derived directly from fits to 
a given set of observations. However, because of their lack of 
theoretical foundation, these formulas cannot be extrapolated to 
explain the properties of galaxies other than those from which 
they are derived. Halfway between these two options is the pos- 
sibility of using analytical expressions endowed with a physical 
basis that enables their application to a wide range of galactic 
and halo parameters. It is precisely with this aim that we here 
introduce a self-consistent pure disk-formation model that fol- 
lows the well-known approach bv iMMWl adapted to the canon- 
ical ACDM concordance cosmology and to four different mass- 
density distributions for the protogalactic dark halos. This up- 
dated scenario is capable of matching simultaneously with very 
good accuracy the zero points and slopes of the observed corre- 
lations in the RVL space of disk galaxies from reasonably real- 
istic values of its input parameters. Yet its most valuable char- 
acteristic is its ease of implementation, as we approximated the 
model predictions for the scale length and rotation speed of disks 
by analytical expressions. The supplied equations can come in 
handy for situations that require the generation of large numbers 
of galaxies with intrinsic attributes in good agreement with the 
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mean observed trends, especially when the relative abundances 
of these objects are known in advance. 

2. Model components 

We recap here the key assumptions and associated equations of 
our self-consistent ACDM-model of disk-galaxy formation: 

1 . In the protogalactic state, the (hot) baryons and dark matter 
are well mixed within virialized spherical halos. Both com- 
ponents have the same distribution of specific angular mo- 
mentum. 

The total angular momentum 7vii of a galactic halo of virial mass 
Mvir is commonly characterized in terms of the dimensionless 
spin parameter 



A^ 



V27;virVv 



-/. 



1/2 



(1) 



which, according to the results of A^-body simulations, follows 
a lognormal distribution with me dian /Ip lying in the range 
0.03 ^ /lo ^ 0.05 (see, for instance. IShaw et al]|2006l and refer- 
ences therein), nearly independent of cosmology, halo environ- 
ment, andredshift (e.g. Lemson & Kauffmann, 1999). In Eq. ([TJ, 
.Rvir is the virial radius inside which the halo mean density, p^j^, 
is Avir times the mean density of the universe at the redshift 
of observation, V^- - GM^ii/Ryi^, and fc is a dimensionless 
function of the halo concentration (see below) that measures 
the deviation of the protogalactic halo's energy from that of a 
sing ular iso thermal sphere with the same mass, -(I /2)MyiiV^. , 
(see iMMWl) . For the family of flat (Q^ h- Qa = 1) cosmogo- 
nies, Avi,(z) ^ {18;r^ -H 82[n(z) - 1] - 39[Q(z) - 1]^}/Q.(z) 
dBrvan & Normanill998l) . 

The halo concentration parameter, c, characterizes the over- 
all shape of a halo density profile by measuring the ratio be- 
tween its outer radius and inner scalelength. Originally intro- 
duced for the Navarro-Frenk- White function, its mean values are 
strongly correlated w ith the halo mass given a cosmology (e.g. 
iNavarro et al.L Il997h . We approximate the mean concentration- 
mass relation at z = in the range of halo masses of interest, 
10'° ^ M/(h'^ Mq) ^ 10", by t he best-fitting power-law rela- 
tion recently inferred bv lMacci6 et al. ( 2008) from relaxed halos 
simulated in the Wilkinson Microwave Anisotropy Probe 5 years 
results (WMAP5) cosmology 



c(Mvir,0) = 9.35 



Mv 



n -0.094 



1012/2-1 Mo 



(2) 



where c = R^ilr-i is defined here in a profile-independent 
form by adopting as the inner characteristic radius of the 
halo density profile the radius r_2, at which its effective log- 
arithmic density slope y = dlnp(r)/dlnr equals -2. The 
WMAPS cosmological parameters are (Qq, JIa.o, ^b,o, h, ct'^, n) - 
(0.26, 0.74, 0.044, 0.72, 0.8, 0.96), implying that p^j^ ~ 96 times 
the critical density for closure at the current epoch. 

2. Disks form smoothly out of cooling flows preserving the spe- 
ciflc angular momentum of the baryons. The cold gas settles 
in centrifugal equilibrium at the center of the halo 's poten- 
tial well following an exponential distribution. 

The fraction of baryons that collect into the central galaxy 
(in the form of stars + cold gas) is defined as 



md = Md/Mv 



(3) 



where the values of this parameter, for which a plausible upper 
limit is the universal baryon fraction /b,o = fib.o/f^o - 0. 17, do 
not s eem to depend much on the halo mass or spin (Sale s et al.L 
|2009). Similarly, the angular momentum of the disk, expressed 
in units of that of its surrounding halo, can be written as 



jd - Jd/Jv 



(4) 



The common yet uncertain assumption that the specific angular 
momenta of the central disk galaxy and of the halo hosting it are 
equal, Jd/Md = 7vii7-Wvii-, is equivalent to setting /'d = m^. 

On the other hand, a thin exponential mass distribution of 
total mass Md, surface density I.(R) - Md/(27r7?2)exp(-/?//?d), 
and a rotation curve V{R), has a total angular momentum 



7d = 2MaRaV, 



1 r 

'2 Jo 



2^-uV(uRd) 



u e 



Vvi 



du = 2Md/?dVv„-/; 



Rd 



(5) 



where the factor /^d is unity for a disk with a flat rotation curve at 
the level Vyir, and where the total circular speed is computed by 
summing in quadrature the contributions from the disk of cold 
baryons, Vd, and from the dark halo, Vh, 



vHR) = Vl(R) ■ 



Vi(r)r- 



(6) 



with R t he cylindrical radius . An e xpression for V^{R) can be 
found in Binney & Tremaine ('2008'), p.lOl, Eq. (2.165). 

Substituting Eqs. ([TJ, ^, and ^ into Eq. ([Sj, one can then 
obtain the disk scalelength as a function of the model parameters 



Rd 



.// 



Rd 



AdRy 



(7) 



with /Id = (jd/'nd)A the effective spin of the disk. 



3. The halo contracts adiabatically and without shell crossing 
to gas inflow. 

According to the adiabatic compression paradigm 

(Blumenthal et al.L Il986h . for a spherical halo in which all 
particles move on circular orbits with velocity V{r), any func- 
tion of the specific angular momentum rV{r) is an adiabatic 
invariant. Assuming that the initial and final mass distributions 
of the different components are spherically symmetric, this 
invariance leads to the relationship 



[Md(r) + M,(r,)(l - md)]r = Mi(n)n , 



(8) 



where r, and r are, respectively, the initial and final radius of 
the spherical shells, Mi(r) is the initial protogalactic halo mass 
profile, and Md(r) comes from the replacement of the final thin 
exponential disk configuration by the spherical density profile 
that has the same enclosed mass. 

The contribution to the total rotation curve (Eq. ||6l) from the 
dark matter (and the remaining hot baryons) is therefore 



y2(r) = GM,(r,)(l-md)/r. 



(9) 



Taking into account that both halo and disk properties are 
directly proportional to their corresponding virial parameters, 
Eq. (|6]l allows one to express the amplitude of the total rota- 
tion curve at a given number of scalelengths and, in particular, 
its peak value, Vmax, in the compact form (cf . iMMW) 



Vn, 



'^ vir/vi 



(10) 



with /vmax a dimensionless factor that, like /sd, depends on the 
adopted halo density law and on the values of parameters A^, c, 
and md. 



Table 1. Halo profiles. 
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Table 2. Model parameters. 



Profile 


p(s) 


c/c. 


Reference 


EIN^i 
HER 
BUR 
NFW 


exp(l-^'/«) 

8/^/(1 + sf 

4/(1 + ^)/(l + ..2) 

4/s/(l + sf 


2 

0.6573 

1 


Einasto & Haud (1989) 
Hemauist(1990) 
Burkert(1995} 
Navarro et al. (1997) 



3. Model predictions 

We now proceed to tune the free parameters of our disk-galaxy 
formation model to match the scaling relations in RVL space 
observed at z ~ 0. For a given halo virial mass, two are the 
free parameters in our modeling: the disk mass fraction, mj, and 
mass-to-light ratio, Tj = M^/L. This latter quantity is needed 
to convert the predicted disk masses into observed luminosities. 
We do not allow the average effective disk spin to vary freely 
however, but use the condition j^ - »Jd to set it equal to three 
representative values of Aq: 0.03, 0.04, and O.Ofl 

We investigated the performance of our model for the four 
functional forms of protogalactic DM halos listed in Table [T] 
They are among the most representative functions used in the 
literature to describe the equilibrium density profiles of halos 
generated in CDM A^-body simulations. All of them are spherical 
density distributions of the form 



p(r) = psp{s) , 



(11) 



where ps = pir^) and r^ are the characteristic density and scale 
radius of the profile respectively, and p(s) is a dimensionless 
function of the dimensionless radius s = r/r^. 

With the aid of the c(Mvii) relation these expressions can be 
reduced to uniparametriqj density laws in which the halo struc- 
ture is fully determined from Mvir. It can be shown that 



1_ 



^-r 



s p(s)ds , 



(12) 



where the characteristic concentration Cs = Ryn/fs is directly 
related to the profile-independent halo concentration parameter 
c defined in Eq. (|2]i (see Table[T]i. 

In order to constrain our model predictions, we consider a 
subset of the SFI-n- sample (ISpringob et al.L 120071) consisting 
of 649 galaxies also included in the co mpilation of ~ 1300 lo - 
cal field and cluster spiral galaxies by 'Court eau et al.l (l2007h . 
The full SFI-i-H- contains measures of intrinsic rotation velocity 
widths reduced to a homogeneous system based on the 21 cm 
spectral line, W, as well as absolute /-band magnitu des for near 
5000 spiral galaxies, while the dataset by Courtea u et al.l pro- 
vides inclination-corrected estimates of disk scalelengths also in 
the /-band (below both observables and model parameters will 
refer to this near- IR band). 

As stated by ! Catinella et al.l (l2007h . for most intermediate 
and bright disks the width of the global H I profile provides a 
more reliable observational estimate of the peak rotation veloc- 
ity than the widths of Ha rotation curves, at least for objects 
not affected by environmental interactions. This is probably be- 
cause the latter are usually evaluated either at a radius where. 



^ We ignore here a possible dependence of this parameter on halo 
mass (e.g. Berta et al., 2008). 

^ The Einasto r''" model has an additional parameter n controlling 
the curvature of the profile. In our modeling this paramete r is kept fixed 
to n = 6, a value representative of galaxy-sized halos jMerrit et al.L 
I2005h . 



Profile 


^d 


nii 


Td/[A(M/L;)o] 


V,„.x/Ki,- 


EIN^ 


0.03 


0.025 


1.20 


1.4 


EINfi 


0.04 


0.050 


1.50 


1.6 


EIN. 


0.05 


0.080 


1.70 


1.8 


HER 


0.03 


0.020 


1.00 


1.4 


HER 


0.04 


0.035 


1.25 


1.5 


HER 


0.05 


0.055 


1.45 


1.7 


BUR 


0.03 


0.030 


1.50 


1.4 


BUR 


0.04 


0.050 


1.60 


1.6 


BUR 


0.05 


0.075 


1.75 


1.7 


NFW 


0.03 


0.030 


1.30 


1.4 


NFW 


0.04 


0.050 


1.50 


1.6 


NFW 


0.05 


0.080 


1.70 


1.8 



on average, they are still rising (e.g., 2.2Rd), or on the asymp- 
totic part of the optical disk. Accordingly, we adopt the approx- 
imation W/2 ^ Vmax, where Vmax is the maximum width of our 
model total speed curve measured within SR^. 



3.1. Scaling laws 

The distribution of /? as a function of V provides the most effec- 
tive way of determining the value of mi - which for bright galax- 
ies represents to a good approximation the stellar mass fraction 
- that best fits the observations for each one of the values of A^ 
under consideration. To allow for a more robust comparison be- 
tween the model predictions and the data, the RV scaling law 
has been recast in the form of the tighter relation between the 
average specific angular momentum of disks computed from the 
fiducial rotation speed of the galaxies, ij = 2/?dVmax, and Vmax- 
In a log-log scale this relationship is expected to follow a straight 
line with a slope near 2 and a zero point that is a sensitive func- 
tion of md. 

In the upper-left panel of Fig. [1] we show the model relations 
that best fit the barycenter of the data cloud for the four halo 
profiles considered and the central value of /Ij (the best values of 
HJd obtained for each one of the three values adopted for /Ij are 
listed in Col. 3 of Table|2]i. It can be seen from this plot that our 
disk models also reproduce the slope of the observed td - Vmax 
scaling law. We note in passing that on the basis of its location in 
this diagram, the angular momentum and disk scale of the Milky 
Way (MW) are unrepresentative of those of a typical spiral (see 
also Hammer et al., 2007). 

With mj fixed and given that the halo concentration is not 
allowed to vary freely, the most sensitive tuning of the other 
free parameter of the model. Yd, is achieved by normalizing 
the model predictions to the observed VL relation. For the latter, 
which is fully independent of surface brightness (Zwaan et aQ, 
1995i ICourteau & RixL Il999l) . we use the calibration of the 
Tully-Fisher (TF) relations hip corrected f rom observational and 
sample biases calculated bv lMasters et al.l (12006) using 807 clus- 
ter galaxies extracted from the SFI-n- catalog, which we rewrite 
in the form 



Mi - 5 log h^ -20.85 - 7.85 [log(2ymax) - 2.5] 



(13) 



to facilitate the comparison with our model predictions. In 
Eq. (fT3] ). M/,0 = 4.11 mag has been adopted to transform model 
luminosities into absolute magnitudes. As in the former case, the 
upper-right panel of Fig. [T]shows that the predicted VL relations 
(again we show only those inferred using the central value of Ad) 
closely match the slope of the empirical estimate. In this case. 
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the best values of Td have been set by minimizing the resid- 
ual between the model predictions and the observed TF rela- 
tionship over the full available range of velocities. As could be 
expected, the agreement between the model predictions and our 
SFI-n— based comparison sample is fairly good too. This panel 
also illustrates the well-known defic iency in luminosity of the 
MW with respect to the TF relation ( Portinari et al.Ll2007h . 

The excellent agreement between predictions and observa- 
tions in the RV and VL planes is maintained for the joint distribu- 
tion of the three variables. The lower-left panel of Fig.[T]depicts, 
again for the central value of A^, the scatter diagram of central 
disk surface density, Eq = M^lilnR^^, and rotation speed. We 
have converted Springob et al. 's data on M/ into total disk lumi- 
nosities, which in turn have been transformed into disk masses 
using the values of Tj derived from the normalization of the VL 
relation. It can be seen that our model predictions are once more 
comfortably close to both the normalization and, in this case, 
curved mean trend delineated by the data. 



3.2. Fitting functions for gaiaxy scaling parameters 

By using the values quoted in Cols. 3 and 4 of Table [2] it is 
straightforward to calculate the average luminosity of a nearby 
disk embedded in a halo of given Myh and A^. However, as 
shown in Sect. |2] each of the remaining fundamental disk prop- 
erties, the characteristic scale and rotation speed, participates 
in the calculation of the other. As a result, they can only be 
computed by applying an iterative procedure that, despite its 
fast convergence, remains cumbersome. For this reason, it is 
very convenient to approximate the dimensionless factors /sj 
and /vmax appearing in the calculation of R^ and Vmax (eqs. Q 
and IfTOl . respectively) by fitting functions. Drawing inspiration 
from MMW, we propose the following fitting formulas, which 
are valid for any of the protogalactic halo mass density profiles 
explored: 

ni+a2»!d/(l+'ld) 



fRd(M, c, nid) ^ ( — — j (fl3 + a^mi + asm\) 



X [flfe + a-jc + a'ijc + fl9/(cmd)] 

1 ,j 0.001 /mj+iimj 



(14) 



/vmaxC^d, C, md) ~ [rrT^] (^2 + ^3'«d + M^d) 

X (^5 +b^c + bilc) . (15) 

The values of the coeflicients corresponding to each profile, 
which are independent of the adopted c(Mvir) relation, are listed 
in Table [3] Both approximations are accurate to within 8 % for 
5 < c < 25, 0.03 < /Id < 0.05, and 0.01 < m^ < 0.1. We note 
that attempts to fit these dimensionless factors using either solely 
polynomials or a linear combination of power laws of parame- 
ters /Id, c, and nii have required a substantially larger number of 
terms to achieve a similarly satisfactory match. Given the count- 
less number of real functions that can be implemented, it would 
be obviously possible to find other formulae also producing good 
fits, but most likely they will be more compUcated than the above 
expressions. 

4. Discussion and conclusions 

We formulated a standard formation model of disk galaxies in- 
side DM halos within the concordant ACDM cosmology that 
simultaneously predicts with high accuracy the main trends of 
the observed fundamental scaling relations of nearby galaxies 
in RVL space. This modeling has been developed with the sole 



aim of deriving physically sound analytical expressions for pre- 
dicting the central properties that characterize the light profiles 
and rotation curves of typical spirals. We supply formulas for 
Einasto r'^^, Hernquist, Burkert, and Navarro-Frenk- White pro- 
togalactic halo mass density distributions that provide a similarly 
good overall description of the data on disks for realistic enough 
values of the model free parameters. We find that, for a given 
/Id, the predictions of the Einasto r'^*, Hernquist, and Navarro- 
Frenk- White profiles are relatively similar, while the Hernquist 
profile - the only density law investigated that does not follow a 
p oc r^^ behavior near 7?vii - requires values of m^ and Td that 
are lower by about a factor of 0.70 and 0.85, respectively. 

The reader may have noticed that our best models yield for 
Td, i.e. for the inverse of the average star-formation efficiency, 
values somewhat lower than those inferred from population syn- 
thesis calculations (e.g. Pizagno et al., 2005). We stress however 
that the observational estimates of this parameter are affected 
by considerable uncertainties, our prediction that the average 
mass-to-light ratio of disks is ~ 1(M /L/)q, which is c onsis- 
tent with submaximal disks arguments ('Courte au & RixL 119991 : 
DCuzio de Naray et al., 2008), as well as relatively close to the 
values a dopted as input in m ore sophisticated models of disk for- 
mation toutton et al.L 120071) . On the other hand, we find that the 
predicted values of md are directly correlated with those adopted 
for /Id. In particular we note that a value of A^ - 0.03, which co- 
incides with the median of the distribution of th e spin parameter 
for relaxed halos derived bv lMaccio et al.l (l2008h . implies a small 
current average galaxy-formation efficiency, md//b,o < 0.2. This 
agrees well with th e predictions of ga laxy evolution from halo 
occupation models feh eng et al.L |2007|) an d methods that match 
the st ellar mass function to that of the halo (IConrov & Wechslen. 
I2009h . Further recent support for low nid (and Ad, according to 
our model) comes fo r insta nce from weak lensing measurements 
(Mandelbaum et al.'^ l2006h and from the roughly universal dis- 
tributions of this parameter obtained by S ales et al.l (l2009h for 
various implementations of feedback in large cosmological A^- 
body/gasdynamical simulations. Notice also the fifth column in 
Table |2] where we list the ratio Vmax/Vvir calculated for a MW- 
mass halo, which increases with increasing A^ and decreasing 
Mvir. As stated bv iDutton et al.l (i2007). the relatively high val- 
ues we predict for this ratio - a characteristic common to stan- 
dard models - would likely hamper a simultaneous match to the 
galaxy LF that, according to semi-analytical models of galaxy 
formation, requires the condition V^ax ~ V'vir- 

We made no attempt to explore the scatter of the observed 
scaling relations and the covariance that exists between model 
parameters, except for Fig. [T] where we carry out a naive com- 
parison between the spread of the data and that resulting from 
taking into account the predicted scale of the probability distri- 
bution of the halo concentration. Including this and other sources 
of scatter, such as the variance of the halo spin parameter, or the 
dependence of the concentration-m ass relation on the adopted 
cosmology (e.g. lMaccio et al.Ll2008l) . would undoubtedly enrich 
the analysis. Yet, a thorough investigation of scatter requires 
dealing with the joint probability distribution of all the param- 
eters entering the model and, in particular, with all their co- 
variances (not just the variances), which ideally should be cor- 
rected from measurement errors. This far exceeds the scope of 
our present research. We note in addition that efforts in the di- 
rection just outlined will soon be much more effective when they 
can be applied to objective, homogeneous, and complete RVL 
datasets free of nontrivial selection biases, as those build from 
the cross-correlation of wide-area spectroscopic optical and H I 
surveys (e.g., Toribio et al. 2010, in preparation). 
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Table 3. Coefficients of the approximations. 



Profile 


fli 


02 


oa 


a4 


as 


as 


ay 


ag 


a. 


EINjs 


-0.041 


5.95 


0.264 


-1.45 


3.35 


3.16 


-0.0311 


5.58 


_ 


HER 


-0.037 


4.86 


0.252 


-1.73 


6.14 


3.13 


-0.0429 


7.62 


- 


BUR 


-0.058 


5.50 


0.275 


-1.20 


0.70 


2.63 


-0.0234 


7.44 


0.027 


NFW 


-0.056 


6.18 


0.267 


-1.46 


3.31 


3.19 


-0.0310 


5.56 


- 




fei 


bi 


b3 


b4 


bs 


bf, 


bi 






EIN. 


-5.50 


0.634 


5.24 


0.78 


1.44 


0.0343 


-0.72 






HER 


-5.58 


0.792 


6.78 


0.16 


1.35 


0.0351 


-1.77 






BUR 


-4.24 


0.722 


4.91 


2.24 


1.32 


0.0394 


-1.03 






NFW 


-5.38 


0.651 


5.38 


0.84 


1.41 


0.0346 


-0.82 







Finally, we wish to comment on the possibility of extending 
our model predictions to distant galaxies by adopting a c(Mvir) 
relationship of the form 



c(Mvir,z) = c(Mvir,0)^(z)- 



(16) 



with c(Mvir, 0) given in Eq. (|2|i and giz) the concentration growth 
factor, which can be taken proportion al to J/(z)^^^ , as fou nd in a 
recent modification of the original | Bullock et alJ (1200 ih model 
for WMAP cosmologies bv lMaccio et all (l2008l) 
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Fig. 1. Scale relations for nearby disks. Upper-left: Disk-specific angular momentum as a func tion of Vmax- The cro ss shows the 
barycenter of the data. Upper-right: /-band TF relation. The dashed line shows the one derived by [Masters et al.l (12006*) from SFI-i-H- 
data (Springob et al., 2007). The cross shows the barycenter of the data. Lower-left: Central-disk surface-mass density vs. Vniax- 
Squares with error bars show the median observational values and the first and third quartiles in each velocity bin. In all these plots 
the solid lines show model predictions for pure disks with A^ = 0.04 using the mean c(Mvir) relation at z = 0, while the dotted curves 
show from top to bottom the predictio ns resulting from a dopting the 2.3th and 97.7th percentiles of the concentration distribution 
(2cr scatter) assuming cringe = 0.11 (iMaccio et al.L 1200 8). The error boxes represent, for comparison, the values measured for 
the MW from: R^ = 2.5 + 0.5 kpc, Md = 5 + 0.5 x 10"' Mq, and V^^^ = 220 + 20 km s"' (Binnev & Tremainei i2008) . and 
L/ = (4 + 0.8) X 10'° Lo (Portina ri et al.L 120071) . The data clouds are build on the compilation o f /-band absolu t e mag nitudes and H I 
rotation widths bv Springob et al] (l2007h and on the /-band disk scalelength measurements bv lCourteau et al.l (l2007l) . 



